Crosslinker energy landscape effects on dynamic mechanical properties of ideal polymer hydrogels.

Reversible crosslinkers can enable several desirable mechanical properties, such as improved toughness and self-healing, when incorporated in polymer networks for bioengineering and structural applications. In this work, we performed coarse-grained molecular dynamics to investigate the effect of the energy landscape of reversible crosslinkers on the dynamic mechanical properties of crosslinked polymer network hydrogels. We report that, for an ideal network, the energy potential of the crosslinker interaction drives the viscosity of the network, where a stronger potential results in a higher viscosity. Additional topographical analyses reveal a mechanistic understanding of the structural rearrangement of the network as it deforms and indicate that as the number of defects increases in the network, the viscosity of the network increases. As an important validation for the relationship between the energy landscape of a crosslinker chemistry and the resulting dynamic mechanical properties of a crosslinked ideal network hydrogel, this work enhances our understanding of deformation mechanisms in polymer networks that cannot easily be revealed by experiment and reveals design ideas that can lead to better performance of the polymer network at the macroscale.

Pathophysiology of chikungunya virus infection associated with fatal outcomes.

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that causes acute, subacute, and chronic human arthritogenic diseases and, in rare instances, can lead to neurological complications and death. Here, we combined epidemiological, virological, histopathological, cytokine, molecular dynamics, metabolomic, proteomic, and genomic analyses to investigate viral and host factors that contribute to chikungunya-associated (CHIK) death. Our results indicate that CHIK deaths are associated with multi-organ infection, central nervous system damage, and elevated serum levels of pro-inflammatory cytokines and chemokines compared with survivors. The histopathologic, metabolite, and proteomic signatures of CHIK deaths reveal hemodynamic disorders and dysregulated immune responses. The CHIKV East-Central-South-African lineage infecting our study population causes both fatal and survival cases. Additionally, CHIKV infection impairs the integrity of the blood-brain barrier, as evidenced by an increase in permeability and altered tight junction protein expression. Overall, our findings improve the understanding of CHIK pathophysiology and the causes of fatal infections.

Thermostability Enhancement of GH 62 α-l-Arabinofuranosidase by Directed Evolution and Rational Design.

GH 62 arabinofuranosidases are known for their excellent specificity for arabinoxylan of agroindustrial residues and their synergism with endoxylanases and other hemicellulases. However, the low thermostability of some GH enzymes hampers potential industrial applications. Protein engineering research highly desires mutations that can enhance thermostability. Therefore, we employed directed evolution using one round of error-prone PCR and site-saturation mutagenesis for thermostability enhancement of GH 62 arabinofuranosidase from Aspergillus fumigatus. Single mutants with enhanced thermostability showed significant ΔΔG changes (<-2.5 kcal/mol) and improvements in perplexity scores from evolutionary scale modeling inverse folding. The best mutant, G205K, increased the melting temperature by 5 °C and the energy of denaturation by 41.3%. We discussed the functional mechanisms for improved stability. Analyzing the adjustments in α-helices, ß-sheets, and loops resulting from point mutations, we have obtained significant knowledge regarding the potential impacts on protein stability, folding, and overall structural integrity.

MassCCS: A High-Performance Collision Cross-Section Software for Large Macromolecular Assemblies.

Ion mobility mass spectrometry (IM-MS) techniques have become highly valued as a tool for structural characterization of biomolecular systems since they yield accurate measurements of the rotationally averaged collision cross-section (CCS) against a buffer gas. Despite its enormous potential, IM-MS data interpretation is often challenging due to the conformational isomerism of metabolites, lipids, proteins, and other biomolecules in the gas phase. Therefore, reliable and fast CCS calculations are needed to help interpret IM-MS data. In this work, we present MassCCS, a parallelized open-source code for computing CCS of molecules ranging from small organic compounds to massive protein assemblies at the trajectory method level of description using atomic and molecular buffer gas particles. The performance of the code is comparable to other available software for small molecules and proteins but is significantly faster for larger macromolecular assemblies. We performed extensive tests regarding accuracy, performance, and scalability with system size and number of CPU cores. MassCCS has proven highly accurate and efficient, with execution times under a few minutes, even for large (84.87 MDa) virus capsid assemblies with very modest computational resources. MassCCS is freely available at https://github.com/cces-cepid/massccs.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study.

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C4-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.

The Role of the Extrafibrillar Volume on the Mechanical Properties of Molecular Models of Mineralized Bone Microfibrils.

Bones are responsible for body support, structure, motion, and several other functions that enable and facilitate life for many different animal species. They exhibit a complex network of distinct physical structures and mechanical properties, which ultimately depend on the fraction of their primary constituents at the molecular scale. However, the relationship between structure and mechanical properties in bones are still not fully understood. Here, we investigate structural and mechanical properties of all-atom bone molecular models composed of type-I collagen, hydroxyapatite (HA), and water by means of fully atomistic molecular dynamics simulations. Our models encompass an extrafibrillar volume (EFV) and consider mineral content in both the EFV and intrafibrillar volume (IFV), consistent with experimental observations. We investigate solvation structures and elastic properties of bone microfibril models with different degrees of mineralization, ranging from highly mineralized to weakly mineralized and nonmineralized models. We find that the local tetrahedral order of water is lost in similar ways in the EFV and IFV regions for all HA containing models, as calcium and phosphate ions are strongly coordinated with water molecules. We also subject our models to tensile loads and analyze the spatial stress distribution over the nanostructure of the material. Our results show that both mineral and water contents accumulate significantly higher stress levels, most notably in the EFV, thus revealing that this region, which has been only recently incorporated in all-atom molecular models, is fundamental for studying the mechanical properties of bones at the nanoscale. Furthermore, our results corroborate the well-established finding that high mineral content makes bone stiffer.

Devising Bone Molecular Models at the Nanoscale: From Usual Mineralized Collagen Fibrils to the First Bone Fibers Including Hydroxyapatite in the Extra-Fibrillar Volume.

At the molecular scale, bone is mainly constituted of type-I collagen, hydroxyapatite, and water. Different fractions of these constituents compose different composite materials that exhibit different mechanical properties at the nanoscale, where the bone is characterized as a fiber, i.e., a bundle of mineralized collagen fibrils surrounded by water and hydroxyapatite in the extra-fibrillar volume. The literature presents only models that resemble mineralized collagen fibrils, including hydroxyapatite in the intra-fibrillar volume only, and lacks a detailed prescription on how to devise such models. Here, we present all-atom bone molecular models at the nanoscale, which, differently from previous bone models, include hydroxyapatite both in the intra-fibrillar volume and in the extra-fibrillar volume, resembling fibers in bones. Our main goal is to provide a detailed prescription on how to devise such models with different fractions of the constituents, and for that reason, we have made step-by-step scripts and files for reproducing these models available. To validate the models, we assessed their elastic properties by performing molecular dynamics simulations that resemble tensile tests, and compared the computed values against the literature (both experimental and computational results). Our results corroborate previous findings, as Young's Modulus values increase with higher fractions of hydroxyapatite, revealing all-atom bone models that include hydroxyapatite in both the intra-fibrillar volume and in the extra-fibrillar volume as a path towards realistic bone modeling at the nanoscale.

Mechanism and biomass association of glucuronoyl esterase: an α/ß hydrolase with potential in biomass conversion.

Glucuronoyl esterases (GEs) are α/ß serine hydrolases and a relatively new addition in the toolbox to reduce the recalcitrance of lignocellulose, the biggest obstacle in cost-effective utilization of this important renewable resource. While biochemical and structural characterization of GEs have progressed greatly recently, there have yet been no mechanistic studies shedding light onto the rate-limiting steps relevant for biomass conversion. The bacterial GE OtCE15A possesses a classical yet distinctive catalytic machinery, with easily identifiable catalytic Ser/His completed by two acidic residues (Glu and Asp) rather than one as in the classical triad, and an Arg side chain participating in the oxyanion hole. By QM/MM calculations, we identified deacylation as the decisive step in catalysis, and quantified the role of Asp, Glu and Arg, showing the latter to be particularly important. The results agree well with experimental and structural data. We further calculated the free-energy barrier of post-catalysis dissociation from a complex natural substrate, suggesting that in industrial settings non-catalytic processes may constitute the rate-limiting step, and pointing to future directions for enzyme engineering in biomass utilization.

Importance of Water in Maintaining Softwood Secondary Cell Wall Nanostructure.

Water is one of the principal constituents by mass of living plant cell walls. However, its role and interactions with secondary cell wall polysaccharides and the impact of dehydration and subsequent rehydration on the molecular architecture are still to be elucidated. This work combines multidimensional solid-state 13C magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) with molecular dynamics modeling to decipher the role of water in the molecular architecture of softwood secondary cell walls. The proximities between all main polymers, their molecular conformations, and interaction energies are compared in never-dried, oven-dried, and rehydrated states. Water is shown to play a critical role at the hemicellulose-cellulose interface. After significant molecular shrinkage caused by dehydration, the original molecular conformation is not fully recovered after rehydration. The changes include xylan becoming more closely and irreversibly associated with cellulose and some mannan becoming more mobile and changing conformation. These irreversible nanostructural changes provide a basis for explaining and improving the properties of wood-based materials.

Nanoscale Wetting of Crystalline Cellulose.

Cellulose possesses considerable potential for a wide range of sustainable applications. Nanocellulose-based material properties are primarily dependent on the structural surface characteristics of its crystalline planes. Experimental measurements of the affinity of crystalline nanocellulose surfaces with water are scarce and challenging to obtain. Therefore, the relative hydrophilicity of different cellulose allomorphs crystalline planes is often inferred from qualitative assessments of their surface and the exposition of polar groups to the solvent. This work investigates the relative hydrophilicity of cellulose surfaces using molecular dynamics simulations. The behavior of a water droplet laid on different crystal planes was used to determine their relative hydrophilicity. The water molecules fully spread onto highly hydrophilic surfaces. However, a water droplet placed on less hydrophilic surfaces equilibrates as an oblate spheroidal cap allowing the measurement of a contact angle. The results indicate that the Iα (010), Iα (11̅0), Iß (010), and Iß (110) faces, as well as the faces of human-made celluloses II and III_I (100), (11̅0), (010), and (110) are all highly hydrophilic. They all have a contact angle value inferior to 11°. Not unexpectedly, the Iα (001) and Iß (100) surfaces are less hydrophilic with contact angles of 48 and 34°, respectively. However, the Iß (11̅0) plane, often referred to as a hydrophilic surface, forms a contact angle of about 32°. The results are rationalized in terms of structure, exposure of hydroxyl groups to the solvent, and degree of cellulose-cellulose versus cellulose-water hydrogen bonds on each face. The simulations also show that the surface oxidation degree tunes the surface hydrophilicity in a nonlinear manner due to cooperative effects involving water-cellulose interactions. Our study helps us to understand how the degree of hydrophilicity of cellulose emerges from specific structural features of each crystalline surface.

Impact of cellulose properties on enzymatic degradation by bacterial GH48 enzymes: Structural and mechanistic insights from processive Bacillus licheniformis Cel48B cellulase.

Processive cellulases are highly efficient molecular engines involved in the cellulose breakdown process. However, the mechanism that processive bacterial enzymes utilize to recruit and retain cellulose strands in the catalytic site remains poorly understood. Here, integrated enzymatic assays, protein crystallography and computational approaches were combined to study the enzymatic properties of the processive BlCel48B cellulase from Bacillus licheniformis. Hydrolytic efficiency, substrate binding affinity, cleavage patterns, and the apparent processivity of bacterial BlCel48B are significantly impacted by the cellulose size and its surface morphology. BlCel48B crystallographic structure was solved with ligands spanning -5 to -2 and +1 to +2 subsites. Statistical coupling analysis and molecular dynamics show that co-evolved residues on active site are critical for stabilizing ligands in the catalytic tunnel. Our results provide mechanistic insights into BlCel48B molecular-level determinants of activity, substrate binding, and processivity on insoluble cellulose, thus shedding light on structure-activity correlations of GH48 family members in general.

QM/MM Simulations of Enzymatic Hydrolysis of Cellulose: Probing the Viability of an Endocyclic Mechanism for an Inverting Cellulase.

Glycoside hydrolases (GH) cleave carbohydrate glycosidic bonds and play pivotal roles in living organisms and in many industrial processes. Unlike acid-catalyzed hydrolysis of carbohydrates in solution, which can occur either via cyclic or acyclic oxocarbenium-like transition states, it is widely accepted that GH-catalyzed hydrolysis proceeds via a general acid mechanism involving a cyclic oxocarbenium-like transition state with protonation of the glycosidic oxygen. The GH45 subfamily C inverting endoglucanase from Phanerochaete chrysosporium (PcCel45A) defies the classical inverting mechanism as its crystal structure conspicuously lacks a general Asp or Glu base residue. Instead, PcCel45A has an Asn residue, a notoriously weak base in solution, as one of its catalytic residues at position 92. Moreover, unlike other inverting GHs, the relative position of the catalytic residues in PcCel45A impairs the proton abstraction from the nucleophilic water that attacks the anomeric carbon, a key step in the classical mechanism. Here, we investigate the viability of an endocyclic mechanism for PcCel45A using hybrid quantum mechanics/molecular mechanics (QM/MM) simulations, with the QM region treated with the self-consistent-charge density-functional tight-binding level of theory. In this mechanism, an acyclic oxocarbenium-like transition state is stabilized leading to the opening of the glucopyranose ring and formation of an unstable acyclic hemiacetal that can be readily decomposed into hydrolysis product. In silico characterization of the Michaelis complex shows that PcCel45A significantly restrains the sugar ring to the 4C1 chair conformation at the -1 subsite of the substrate binding cleft, in contrast to the classical exocyclic mechanism in which ring puckering is critical. We also show that PcCel45A provides an environment where the catalytic Asn92 residue in its standard amide form participates in a cooperative hydrogen bond network resulting in its increased nucleophilicity due to an increased negative charge on the oxygen atom. Our results for PcCel45A suggest that carbohydrate hydrolysis catalyzed by GHs may take an alternative route from the classical mechanism.

Transition Path Sampling Study of the Feruloyl Esterase Mechanism.

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

Insights into the dual cleavage activity of the GH16 laminarinase enzyme class on ß-1,3 and ß-1,4 glycosidic bonds.

Glycoside hydrolases (GHs) are involved in the degradation of a wide diversity of carbohydrates and present several biotechnological applications. Many GH families are composed of enzymes with a single well-defined specificity. In contrast, enzymes from the GH16 family can act on a range of different polysaccharides, including ß-glucans and galactans. SCLam, a GH16 member derived from a soil metagenome, an endo-ß-1,3(4)-glucanase (EC 3.2.1.6), can cleave both ß-1,3 and ß-1,4 glycosidic bonds in glucans, such as laminarin, barley ß-glucan, and cello-oligosaccharides. A similar cleavage pattern was previously reported for other GH16 family members. However, the molecular mechanisms for this dual cleavage activity on (1,3)- and (1,4)-ß-D-glycosidic bonds by laminarinases have not been elucidated. In this sense, we determined the X-ray structure of a presumably inactive form of SCLam cocrystallized with different oligosaccharides. The solved structures revealed general bound products that are formed owing to residual activities of hydrolysis and transglycosylation. Biochemical and biophysical analyses and molecular dynamics simulations help to rationalize differences in activity toward different substrates. Our results depicted a bulky aromatic residue near the catalytic site critical to select the preferable configuration of glycosidic bonds in the binding cleft. Altogether, these data contribute to understanding the structural basis of recognition and hydrolysis of ß-1,3 and ß-1,4 glycosidic linkages of the laminarinase enzyme class, which is valuable for future studies on the GH16 family members and applications related to biomass conversion into feedstocks and bioproducts.

SARS-CoV-2 receptor is co-expressed with elements of the kinin-kallikrein, renin-angiotensin and coagulation systems in alveolar cells.

SARS-CoV-2, the pathogenic agent of COVID-19, employs angiotensin converting enzyme-2 (ACE2) as its cell entry receptor. Clinical data reveal that in severe COVID-19, SARS-CoV-2 infects the lung, leading to a frequently lethal triad of respiratory insufficiency, acute cardiovascular failure, and coagulopathy. Physiologically, ACE2 plays a role in the regulation of three systems that could potentially be involved in the pathogenesis of severe COVID-19: the kinin-kallikrein system, resulting in acute lung inflammatory edema; the renin-angiotensin system, promoting cardiovascular instability; and the coagulation system, leading to thromboembolism. Here we assembled a healthy human lung cell atlas meta-analysis with ~ 130,000 public single-cell transcriptomes and show that key elements of the bradykinin, angiotensin and coagulation systems are co-expressed with ACE2 in alveolar cells and associated with their differentiation dynamics, which could explain how changes in ACE2 promoted by SARS-CoV-2 cell entry result in the development of the three most severe clinical components of COVID-19.

Publisher Correction: Structural insights into ß-1,3-glucan cleavage by a glycoside hydrolase family.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Molecular Modeling of Aquaporins from Leishmania major.

Aquaporins are membrane proteins responsible for permeating water, ions, dissolved gases, and other small molecular weight compounds through the protective cell membranes of living organisms. These proteins have been gaining increased importance as targets for treating a variety of parasitic diseases, since they control key physiological processes in the life cycle of parasitic protozoans, such as the uptake of nutrients, release of metabolites, and alleviation of osmotic stress. In this work, we use homology modeling to build three-dimensional structures for the four main aquaporins encoded and expressed by Leishmania major, a protozoan that causes leishmaniasis and affects millions of people worldwide. Physico-chemical properties of the proposed models for LmAQP1, LmAQPα, LmAQPß, and LmAQPγ are then investigated using molecular dynamics simulations and the reference interaction site model (RISM) molecular theory of solvation. Pore characteristics, water permeation, and potential of mean force across the AQP channels for water, methanol, urea, ammonia, and carbon dioxide are examined and compared with results obtained for a protozoan (Plasmodium falciparum) aquaporin for which a crystal structure is available.

Adsorption of CO2 and CH4 in MIL-47 investigated by the 3D-RISM molecular theory of solvation.

Metal-organic frameworks (MOFs) comprise a class of highly porous nanomaterials formed by the assembly of organic molecular templates connected by metal ions. These materials exhibit a large diversity of pore size and geometry, topology, surface area, and chemical functionality. MOFs are particularly promising materials for developing new technologies for capture and storage of greenhouse gases such as methane and carbon dioxide. Here we apply the three dimensional reference interaction site model (3D-RISM) molecular theory of solvation to study the interactions of CO2 and CH4 with the metal-organic material MIL-47. The 3D-RISM integral equations were solved to determine the three dimensional density correlation functions of the gas (solvent) relative to the atomic positions of the MIL-47 framework, treated as static solute sites. The computed solvent spatial distributions inside MIL-47 pores were used to identify whether or not there exist preferable binding sites and the binding free energy landscape for the gas of interest at low computational costs compared with other molecular modeling techniques, such as grand canonical Monte Carlo and molecular dynamics simulations. The 3D-RISM formalism was applied to pure CO2, pure CH4, and binary mixtures of these gases of various compositions under different pressure conditions. The results indicate that both gases bind very weakly to MIL-47 and that this material exhibits nearly vanishing CO2/CH4 selectivity. The 3D-RISM computations presented here can be extended to investigate the physical adsorption of gases on other MOFs and nanoporous materials, providing an alternative low-cost computational approach to study gas capture and storage in nanoporous materials in general and, in particular, to determine the binding free-energy landscape in these systems.

Structural insights into ß-1,3-glucan cleavage by a glycoside hydrolase family.

The fundamental and assorted roles of ß-1,3-glucans in nature are underpinned on diverse chemistry and molecular structures, demanding sophisticated and intricate enzymatic systems for their processing. In this work, the selectivity and modes of action of a glycoside hydrolase family active on ß-1,3-glucans were systematically investigated combining sequence similarity network, phylogeny, X-ray crystallography, enzyme kinetics, mutagenesis and molecular dynamics. This family exhibits a minimalist and versatile (α/ß)-barrel scaffold, which can harbor distinguishing exo or endo modes of action, including an ancillary-binding site for the anchoring of triple-helical ß-1,3-glucans. The substrate binding occurs via a hydrophobic knuckle complementary to the canonical curved conformation of ß-1,3-glucans or through a substrate conformational change imposed by the active-site topology of some fungal enzymes. Together, these findings expand our understanding of the enzymatic arsenal of bacteria and fungi for the breakdown and modification of ß-1,3-glucans, which can be exploited for biotechnological applications.

Collision Cross Section Calculations Using HPCCS.

A technical overview of the High Performance Collision Cross Section (HPCCS) software for accurate and efficient calculations of collision cross sections for molecular ions ranging from small organic molecules to large protein complexes is presented. The program uses helium or nitrogen as buffer gas with considerable gains in computer time compared to publicly available codes under the Trajectory Method approximation. HPCCS is freely available under the Academic Use License at https://github.com/cepid-cces/hpccs .